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Toward  the  end  of  the  research  period,  exploration  of  a  new  method 
commenced.  This  procedure,  which  appears  to  furnish  the  desired  math¬ 
ematical  simplicity,  will  be  studied  further  during  the  forthcoming 
research  period. 

(d)  Finite  Rotations.  Some  problems,  require  the  proper 
description  of  geometric  nonlinearities.  One  such  problem  concerns 
dowel  action.  In  a  narrow  zone  adjacent  to  the  crack/slip  plane,  the 
dowels  that  penetrate  this  plane  undergo  finite  rotations.  In  order 
to  accommodate  this  effect,  finite  rotations  were  included  in  the 
mixture  formulation.  This  phenomenon  was  modeled  using  a  von  Karman- 
type  approximation  that  is  valid  for  finite  but  moderate  rotations. 
Details  concerning  this  effort  can  be  found  in  paper  No.  12  of 
Section  5. 

4.4.2  Model  Validations 

(a)  Behavior  in  Shear  (the  Dowel  Problem).  A  considerable 
effort  was  made  during  the  research  period  to  validate  the  mixture 
formulation  in  shear.  For  this  purpose  the  dowel  problem,  which  was 
originally  investigated  during  an  earlier  research  period,  was 
examined  again  in  detail.  The  results  of  this  examination  led  to  the 
reformulation  of  the  mixture  model  to  accommodate  finite  rotations  as 
was  noted  above.  Validations  were  then  re-conducted.  These  valida¬ 
tions  consisted  of  comparisons  of  simulated  versus  measured  behavior 
for  the  problem  described  in  subsection  4.4.2.  Details  concerning 
this  investigation  are  contained  in  paper  No.  12  of  Section  5. 
Typical  theoreti ca I -exper imenta I  comparisons  were  presented  in  subsec¬ 
tion  4.4.2.  These  comparisons  indicate  that  the  current  mixture 
description  furnishes  excel  I ent  global  response  simulation  capability 
in  the  above  shear  mode. 


9 


% 


9 


9 


* 


(b)  Behavior  in  Bending.  An  attempt  was  made  to  predict  the 
stiffness  degradation  and  crack  distribution  evaluations  associated 
with  f lexura I -type  loading  of  reinforced  concrete  beams.  This  was 
approached  by  first  deriving  a  "mixture"  beam  theory  from  the  full  3D 
mixture  model  via  suitable  constraints  on  the  global  displacement 
field.  This  is  analogous  to  the  derivation  of  beam  models,  such  as 
the  Bernou II i -Euler  Theory,  from  3D  elasticity. 

The  advantage  of  a  beam-type  model  over  a  full  2D  or  3D  analysis 
is  supposed  to  be  the  ease  with  which  one  can  effect  analytical  or 
numerical  solutions.  However,  serious  complications  were  encountered 
in  this  case  which  led  to  numerical  problems.  The  trouble  was  traced 
to  the  method  by  which  "non-thru"  cracks  were  being  incorporated.  A 
new  procedure  has  been  proposed  and,  during  the  forthcoming  research 
period,  another  attempt  will  be  made  to  examine  the  beam  bending 
example. 
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APPENDIX  A 


NEW  FAILURE  SURFACE  FOR  THE  ENDOCHRONIC  CONCRETE  MODEL 

This  appendix  describes  a  new  failure  surface  for  concrete  and 
the  procedure  for  introducing  it  into  the  endochronic  concrete  model. 
Unlike  the  earlier  failure  surface  used  with  this  model,  the  new 
failure  representation  accounts  for  the  influence  of  the  intermediate 
principal  stress  on  failure,  which  is  known  to  be  important  for  con¬ 
crete,  and  provides  an  excellent  correlation  of  various  concrete 
triaxial  failure  data. 


For  the  purposes  of  the  present  discussion,  let  us  recall  the 
expression  for  deviatcric  stress,  s,  as  given  by  the  model  i.e., 


(A- 1 ) 


where  the  intrinsic  time  for  shear, 


’s’ 


is  defined  by  the  expression 


d  z  = 
s 


dz 

F 

r 


(A— 2) 


Here,  F  the  shear  hardening  function,  is  responsible  for  the 
characteristics  of  the  failure  surface. 

In  the  original  endochronic  concrete  model  (Valanis  and  Read, 
1985,1986),  Fs  was  taken  to  depend  on  the  pressure,  0,  in  a  linear 
manner  This  resulted  in  a  Vohr-Coulomb  type  failure  surface  which 
was  independent  of  the  intermediate  principal  stress  and  which  cut  the 
T-plane  in  a  circle.  Despite  its  simplicity,  this  simple  failure 
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criterion  was  evidently  adequate  for  describing  the  data  from  Sca- 
vuzzo,  et  aL.,  (1983),  which  were  obtained  for  stress  states  that  were 
sufficiently  far  from  failure  that  no  significant  cracking  occurred. 
Under  these  circumstances,  it  is  to  be  expected  that  the  test  results 
will  not  be  sensitive  to  the  details  of  the  failure  surface,  and  this, 
in  fact,  appeared  to  be  the  case. 


To  obtain  accurate  descriptions  of  response  near  failure,  a  more 
sophisticated  failure  criterion  is  required,  specifically  one  that 
reflects  the  influence  of  the  intermediate  principal  stress.  For 
isotropic  media,  the  failure  criterion  in  its  most  general  form,  is 
expressible  in  terms  of  three  independent  stress  invariants.  For  a 
convenient  representation  in  principal  stress  space,  it  is  useful  to 
consider  the  three  stress  invariants  a,  J2  and  6  where 


a  =  tr  (a) 


s :  s 

M  M 


e 


-1 

cos 


3  f3  J3 
2(J2)3/2 


(A-3) 

(A-4) 

(A— 5) 


and  Jg  is  defined  as 
J3  =  3  det  & 


(A— 6) 


The  general  form  of  the  failure  criterion  for  isotropic  media  can 
therefore  be  written  in  terms  of  these  stress  invariants  as 


To  determine  the  relation  between  the  functions  Fs  and  con¬ 
sider  the  case  of  radial  deformation  at  constant  hydrostatic  pressure. 
For  this  case,  we  can  write 


de^  =  n  I  I  de*3  II  =  n  d<T 

a*  n,  «v  A*  j  e 


(A— 8) 


so  that 


^(Zs-Z']d^d2' 


=  C  p(z&-  z')jp  dz- 


(A— 9) 


During  the  initial  hydrostatic  loading  process  for  0  <  zs  <z°s,  we 
have  =  0.  Thus,  if  we  set 


w  =  z  -  z 
s  s 


(A- 10) 


it  fol lows  that 


z  -  z  =  w  ♦  z  -  z  =  w  -  w 
s  s  s  s 


(A— 11) 


As  a  result,  Eq.  (A-9)  can  be  rewritten  in  the  form 


«  '(v  z  )d rdz' 

J0 


tt  dz'  =  jq  -  w')^r  dw' 


(A-12) 


It  has  been  shown  by  Read  (1985)  that  when  the  hydrostatic  hardening 
function  F^  is  of  the  form: 


fh  *  e 


(4-13) 


v'«*  v 


Afe. 


where  p  is  a  positive  constant,  it  follows  that 


dc  F  (2cw  +  c2w2) 
^s  _  _s _ 

dw  (1  +  cw) 


where 


c  =  £  F 
k  s 


(A— 14) 


(A— 15) 


and  p,  k  are  material  constants. 
(A-12)  gives 


Substitution  of  Eq.  (A-14)  into  Eq. 


r 2  cw'  +  c2(w')2 
\  L _ 


S  =  n  F  p(  W  -  w') - 7Z - TT- 

*  ~s  ’  (1+cw) 


(A-16) 


Using  a  result  given  by  Valanis  and  Read  (1984),  it  can  be  shown  that 


[2  cw'  +  c2(w')2] 
(1  cw') 


lim  />{ w  -  w  )  ;  J  dw'  =  Mo 

r  '  fl  ♦  cw  )  « 


(A— 1 7) 


where 


=  Pi*)** 


(A— 18) 


Thus,  from  Eqs.  (A-16)  and  (A-17)  ,  we  can  write 
s^  =  I im  s  =  n  F  M 

**  ^  N  C  OO 

w+* 


(A— 19) 


where  the  superscript  f  denotes  failure.  It  then  follows  that  at 
fail ure  we  have 


The  relation  between  the  functions  Fg  and  ^  is  then  obtained  by  com¬ 
paring  Eqs.  (A-7)  and  (A-20) ,  with  the  result: 


(A— 21) 


Thus,  to  define  Fg,  the  function  4>(o,8)  needs  to  be  specified. 

There  are  a  number  of  advanced  failure  criteria  for  plain  con¬ 
crete  available  which  have  the  general  form  of  Eq.  (A-7)  (Mills  and 
Zimmerman,  1970;  Wiliam  and  Warnke,  1974;  Ottosen,  1977;  Lade,  1982; 
Peyton,  1983;  Podgorsk i ,  1985).  Considering  the  large  scatter  in  data 
between  experimenters  as  well  as  between  different  devices  (Hegemier 
and  Read,  1986),  most  of  these  advanced  failure  criteria  provide  a 
reasonably  accurate  description  of  existing  concrete  failure  data. 
The  reader  is  referred  to  Ottosen  (1977)  and  Podgorsk i  (1985)  for  a 
critical  assessment  of  most  of  these  failure  models. 

Inasmuch  as  most  of  the  above  failure  models  are  reasonably 
accurate,  and  none  appear  to  have  a  significant  advantage  in 
predictive  capability  or  numerical  implementation  over  the  others,  we 
adopted  for  use  in  the  present  study  the  failure  criterion  that  was 
developed  by  S-CUBED  (Peyton,  1983)  during  the  previous  AFOSR  program 
(see  Hegemier,  et  aL,  1983).  A  summary  of  this  criterion  is  given 
below. 


Ift! 


£ 


K** 


The  failure  criterion  of  Peyton  (1983)  is  of  the  form 


(A-22) 


where  r  and  a  are,  in  general,  smooth  monotonic  functions  of  pressure 
a ,  i  .  e . , 


r  =  r(ff) 
a  =  a  (a) 


(A-23) 


The  role  of  the  function  r(a)  is  to  describe  the  effect  of  pressure  on 
the  meridians  of  the  failure  surface,  while  the  function  a(cr)  deter¬ 
mines  the  manner  in  which  the  trace  of  the  failure  surface  in  the  t- 
plane  changes  with  pressure.  In  order  for  the  trace  to  be  convex,  it 
is  necessary  that  0  <  o  <  1.  When  a  is  a  continuous  decreasing  func¬ 
tion  of  pressure,  the  trace  of  the  failure  surface  in  the  ir-plane 
changes  smoothly  from  a  triangle  (a  =  1)  to  a  c i rc I e  (a  =  0)  w i th 
increasing  pressure. 


In  order  to  express  Eq.  (A-22)  in  the  form  of  Eq.  (A-7) ,  we 
introduce  the  stress  invariant  J,  where 


_  iD*  J3 
=  2  (J2)3/2 


(A-24) 


Upon  defining  an  angle  d  in  the  ff-plane  as  shown  in  Figure  A.l,  it 
fol lows  that 

J  =  -  sin  36  (A-25) 
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A 


i f  we  set 


X 


sin  36  =  aJ 


(A— 32) 


Thus,  from  Eq.  (A-32b) ,  we  can  write, 


=  |  [sin  ^  (aJ)  ♦  2nx] 


(A- 33) 


and,  returning  to  the  original  variables,  it  follows  that 


J3 


2  sin 


'l  -  -1 

,n  [3  Sl" 


,  i\  2* 

(aJ)  ♦  -3 


(A-34) 


which  has  the  form  of  Eq.  (A-7) .  To  fit  this  model  to  a  particular 
set  of  data,  the  functions  t(o)  and  a(o)  need  to  be  specified.  A 
procedure  for  doing  this  is  described  below. 

Generally,  failure  data  are  available  from  triaxial  compression 
and  triaxial  extension  tests.  Using  such  data,  we  define  a  ratio  r  as 


r 


compress i on 
extens i on 


(A-35) 


where  both  values  of  D2  are  taken  at  the  same  pressure.  Using  Eq . 
(A-34)  with  J  =  +1  for  compression  and  J  =  -1  for  extension,  it  is 
straightforward  to  show  that  r  has  the  form: 
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This  expression  can  be  solved  for  a  to  give 


a 


sin  1 3  tan  ^13  ^  ^  jj 


(A- 37) 


Therefore,  by  knowing  how  (^2) compress i on  and  (^2) extens i on  depend 
upon  a  from  data,  the  dependence  of  3  on  n  can  be  determined  from  Eq . 
(A-37)  . 


Having  determined  a(a),  and  with  the  variation  of  (■J'^compr 
with  a  known  from  data,  the  function  t(o)  can  be  determined  from  the 
fol lowing  equation 


'W  *  j|  sin(!  sin"X  “  *  ^tH^ccpr.  (A~38) 

which  is  obtained  by  simply  solving  Eq .  (A-34)  for  r(a)  and  setting 

J  =  *1. 

The  above  procedure  was  used  to  determine  the  forms  of  a  and  r 
for  plain  concrete  from  the  failure  data  shown  earlier  in  Figure  3.2. 
Here,  failure  data  from  both  triaxial  compression  and  triaxial 
extension  tests  on  plain  concrete  conducted  by  several  investigators 
on  samples  of  several  different  strengths  are  depicted.^  After 


4.  Note  that  the  data  of  Scavuzzo,  et  ai.  (1983)  increasingly  deviate 
from  the  other  data  shown  in  the  figures  with  increasing  pressure. 


applying  the  above  methods  for  determining  a  and  r  from  the  data,  it 
was  found  that 


a  -  a 


r  =  t  *  B  ia  -  7 
o  ro  'o 


(A-39) 


where 


a  =0.83 
o 


r  =  -  1.913  f' 
o  c 


(A-40) 


p  =  2.i44rr 
ro  c 


7  =  -  0.896  f' 

'o  c 


The  excellent  ability  of  the  model,  with  only  four  material  para¬ 
meters,  to  correlate  the  data  was  shown  earlier  in  Figure  3.2. 


Returning  to  the  endochronic  model  for  plain  concrete  (Valanis 
and  Read,  1985;1986),  we  note  that  the  kernel  function  />( z)  was  taken 
in  the  form: 


Z-a  3 

V  r 


(A— 4 1 ) 


so  that,  in  view  of  Eq.  (A-18)  we  have 


M  =  Z  ~ 

oo  ^ — •  a 

r 


(A-42) 


Upon  substituting  the  values  for  A2  and  ar  given  by  Valanis  and  Read 
(1985; 1986)  into  Eq.  (A-18),  it  follows  that 


•'.V.  -’.V  V  ‘ •' 
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APPENDIX  B 

NEW  DAMAGE-CRACKING  MODEL  FOR  CONCRETE 

The  purpose  of  this  appendix  is  to  describe,  in  detail,  a 
damage-cracking  model  that  is  being  developed  under  the  present 

program  for  ultimate  use  in  conjunction  with  the  endochronic  concrete 
model.  The  model  was  formulated  by  Valanis  (1935)  and  is  in  the  early 
stages  of  development.  The  discussion  given  below  largely  follows 
that  given  by  Valanis  (1985)  and  does  not  reflect  various  topics  that 
are  currently  under  study,  which  include  the  addition  of 

(a)  plasticity,  (b)  non-d i rect i ona I  damage  due  to  volume  compaction, 
(c)  size  effects  and  (d)  strain  rate  dependence. 

Formulation  of  Model 


Consider  a  linearly  elastic  material  which  is  isotropic  in  its 
virgin  unstrained  state  and  undergoes  small  deformation  at  isothermal 
conditions.  As  the  deformation  proceeds,  damage  gradually  develops, 
which  in  turn  reduces  the  material  integrity.  For  sufficiently  large 
deformation,  the  accumulated  damage  can  lead  to  complete  fracture,  in 
which  the  material  cannot  support  tensile  stress  in  a  direction  normal 
to  the  crack.  We  shall  formulate  this  damage-cracking  process  within 
the  context  of  irreversible  thermodynamics  and,  for  this  purpose,  we 
introduce  the  free  energy  per  unit  volume.^  Since  fracture  leads  to 
a  reduction  in  material  integrity,  we  introduce  an  integrity  tensor  ^ 
which  will  be  discussed  at  length  later  and  which  is  symmetrical  such 
that  ^  =  i5  (the  unit  tensor),  when  the  material  is  in  its  virgin 


5.  This  nomenclature  is  preferred  to  free  energy  density,  in  view  of 
the  subsequent  concept  of  "irreducible"  material  volume  More  will 
be  said  about  this  later  on. 
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I  (undamaged)  state,  and  Q  when  the  material  has  fully  failed,  i.e., 

;  when  the  material  cannot  support  stress  in  any  direction.  We  thus 

|  set : 


l 

*  =  *(£,*)  (B-l) 

where  £  denotes  the  strain  tensor.  In  a  thermodynami c  sense,  $  new 
plays  the  role  of  an  internal  variable,  in  which  case  we  can  write  the 
fol lowing  relations: 


a 


(B-2) 


where  a  is  the  stress  tensor  and  Q  denotes  the  internal  force,  dual  to 
$  which  is  driving  the  fracturing  process. 


In  view  of  the  definition  of  the  following  relations  must 

hold: 


*(£,Q)  =  0  ,  *(0,*)  =  o  ,  (B-3) 

the  first  meaning  that  a  fully  failed  material  cannot  contain  free 
energy  and  the  second  that  an  unstrained  material  must  have  zero  free 
energy;  both  conditions  are  relative  to  the  reference  state.  Further¬ 
more,  the  following  conditions  must  also  hold: 

2(e,9)  =  0  ,  g(Q,  £)  =  0 

9 (£»S)  =  9 

Equations  (B-4)  require  that  the  stress  vanish  in  the  fully  failed 
material  and  that,  since  the  material  remains  elastic  in  the  damaged 
state,  the  stress  will  vanish  at  zero  strain.  Equation  (3.5)  stipu¬ 
lates  that  the  internal  fracture  causing  force  Q  must  vanish  in  the 
fully  failed  material. 


(B-4) 

(B-5) 


2 


If  we  now  expand  41  in  a  Taylor  senes  in  £,  retain  terms  no 
higher  than  the  quadratic,  and  observe  relations  (3  3),  we  obtain  the 
f  O  I  lowing  form  of  : 


i  j  £  k£. 


(B-6) 


where  C  =  C(^) ,  which  by  virtue  of  Eqs .  (B- 3)  may  contain  no  constant 
terms  nor  linear  terms  in  $ .  Since  the  material  is  assumed  to  be 
isotropic  in  its  virgin  state,  C  way  be  represented  in  terms  of  outer 
products  of  the  unit  tensor  £  and  Furthermore,  since  C  is  purely 
quadratic  in  0,  it  has  no  other  representation  than 

cijkt*xhj^kt*2',^kV  <B-7> 


in  view  of  the  symmetries  imposed  upon  it  by  the  symmetry  of  fc  implied 
in  Eq.  (B-6).  The  constants  X  and  2)i  in  Eq.  (B- 7)  must  indeed  by  the 
Lame  constants  of  the  virgin  material,  since 

C  i  j  kC.  =  X  6\j6kl  *  2fi6,k6}H  (B"8) 
when  ^  =  (5 .  Thus,  we  can  write 


*  *  I  Mh/ij)2  *  r  *,kV,;E 


k  l  ■ 


(B-9) 


so  that,  in  view  of  Eq .  (B-2a)  ,  a.  is  given  by  the  expressii 


aij  =  X  t]jhlLekll  *  tk^l^kl 


(B- 1 0) 


This  equation  defines  the  const1 tuti ve  response  of  an  elastic  frac¬ 
turing  material,  once  we  establish  how  the  integrity  tensor  ^  is 


•At. 


related  to  the  fracturing  process.  In  particular,  an  equation  is 
needed  which  describes  the  manner  in  which  ^  evolves  with  accumulated 
damage,  and  this  is  considered  below. 

Evolution  Equation  for  & 

When  the  materia!  is  not  strain  rate  sensitive,  an  increment  of 
tensile  strain  will  produce  an  increment  of  damage.  Specifically, 
given  an  increment  in  strain,  d£,  let  dea  be  its  eigenvalues  and  n;fl 
its  eigenvectors.  If  d£  and  £  are  coaxial,  then  dea  will  constitute 
an  increase  in  tensile  strain  if  dea  >  0  and  e0)  0,  where  ea  denote 
the  eigenvalues  of  £.  In  the  more  general  case  where  d£  and  £  are  not 
co-axial,  dea  is  said  to  be  an  increase  in  tensile  strain  if  dea  >  0 
and  £jjnjanja  =  efln  >  0  (fl  not  summed).  Now  let  d$  be  the  change  in 
the  integrity  tensor  ^  due  to  the  increment  d£,  and  let  be  the 
eigenvector  of  d$j.  We  now  specify  that  d $  and  d£  are  co-axial,  so 
that 


a  a 
n,  =  n 


(B-l 1) 


The  physical  significance  of  Eq .  (B-ll)  is  examined  in  the  discussion 
following  Eq.  (B-14) . 

To  completely  characterize  the  evolution  of  the  relation 
between  d^a  and  dea  needs  to  be  specified.  For  this  purpose,  we  adopt 
the  following  expression: 


kde  ,  for  de  >  0  and  £  >  0 

n  — 

0  ,  otherwise 


(B- 13) 


Here,  k  and  m  are  positive  constants,  and  ^  i  j n? n<2  j  (a  not 

summed) . 

In  view  of  Eq .  (B-ll),  we  can  therefore  write 


, 

,a  a  a 
\9  n. n .  d{ 
lrnJ  i  j  sc 


(B- 14) 


With  this  evolution  equation  for  l,  the  constitutive  description  of 
the  elastic  fracturing  solid  is  now  complete,  once  the  material  para¬ 
meters  X,  n,  k  and  m  have  been  specified. 

Physical  Interpretation  of  Equations 

The  physical  meaning  of  the  above  equations,  and  their  relation 
to  the  fracturing  process,  is  best  understood  by  reference  to  Eq .  (B- 
10),  which  relates  the  stress  tensor  to  the  strain  tensor.  Let  be 
the  eigenvectors  of  ^  and  let  a  and  £  be  referred  to  a  system  of 
coordinates  Then,  the  components  of  a  and  £  on  the  XQ  coordinate 
system  can  be  written  as: 


=  a- 


'  J 


•  NrNs 


rs 


=  e . 


.NrNs 
•J  i  J 


(B— 15) 


(B-16) 


In  this  system  of  coordinates,  ^  is  diagonal,  i.e.,  ^rs  =  0  for  r  t  s, 
and  Eq.  (B-10)  takes  the  form: 


o  =  X  6  4  e  +  2u  4  4  e  (B-l 7) 

rs  rrsrmn  mn  r  rrmrsn  mn  v  ' 

Since  £  is  diagonal,  we  can  write: 

ffll  =  X  ^llX^r£rr  +  *llell 

r=l 

°21  =  °12  =  2m  ^11^22C12  (B-18) 


^31  =  ^13  =  2/i  ^11^22e13 

Let  us  now  single  out  the  eigenvector  N*  and  examine  what  happens  to 
the  stress  when  the  eigenvalue  ^  (=  ^j)  vanishes.  Note  from  Eqs .  (B- 
18)  that,  if  =  0,  it  follows  that  =  a^2  =  ^13  =  0#  that  is, 
the  solid  cannot  support  stresses  on  a  plane  normal  to  N*,  on  which 
^  =  0.  The  physical  meaning  of  this  is  that  the  decrease  of  ^ 
represents  damage  norma  I  to  the  plane  of  N^  and,  as  such,  it  is  a 
measure  of  the  plane  microcracks  that  have  developed  normal  to  N* 
and/or  the  increase  in  the  size  of  such  cracks,  so  that  when  =  0, 
a  plane  crack  in  the  accepted  sense  has  formed  across  the  material 
element  on  a  plane  with  normal  N^,  so  that  the  element  cannot  support 
tensile  or  shear  stresses  on  that  plane. 

On  the  basis  of  the  above  observations,  the  physical  meaning  of 
Eqs.  (B-l 1 ) ,  (B-12)  and  (B-13)  now  becomes  clear.  An  increase  in 
tensile  strain  de°  on  a  plane  with  normal  na  causes  planar  damage  d $a 
on  that  plane,  where  d^a  is  given  by  Eq.  (B-12),  in  accordance  with 


the  observation  that  "planar  microcracks  form  perpendicular  to  the 
direction  of  the  principal  tensile  strain." 

The  constant  k  in  Eq.  (B-13)  is  a  material  parameter  that 
reflects  the  fracture  resistance  or  fracture  toughness  of  the  mate¬ 
rial.  Eq.  (B-12)  is  representative  of  processes  involving  annihila¬ 
tion  of  populations  and  catastrophic  processes  in  systems  where  the 
increment  of  annihilation  is  proportional  to  the  state  of  integrity, 
which  is  given  here  by  the  factor  This  form  is  representative 
of  such  systems,  even  though  any  monoton i ca I  I y  decreasing  function 
f ,  such  that  f (0)  =  0,  will  likely  do  as  we  I  I .  Finally,  we  note 
that  the  damage  tensor  Q  can  be  expressed  in  terms  of  ^  as  follows: 


D..  =  H 
'•>  a 


nV 


(B— 19) 


The  physical  interpretation  of  the  governing  equations  is  now 
complete. 

3.3.3  Application  to  Some  Simple  Cases  of  Homogeneous  Deformation 

In  this  section,  analytic  solutions,  based  upon  the  constitutive 
model  described  above,  are  presented  for  several  simple  cases  of  homo¬ 
geneous  deformation  including  simple  tension,  simple  compression  and 
simple  shear.  These  solutions  provide  valuable  insight  into  the  char¬ 
acteristics  of  the  model  and  reveal  its  remarkable  predictive  capabil¬ 
ities,  despite  the  fact  that  it  involves  only  four  material 
parameters . 


As  a  preliminary  development,  consider  the  case  in  which  the 
strain  field  is  homogeneous  and  consists  only  of  principal  strains, 
i.e.,  e jj  =0  for  i  i  j.  Thus,  d£  and  £  are  always  coaxial  and  $  is 
coaxial  with  £  by  virtue  of  Eq.  (B-ll) .  As  a  result,  ^  is  diagonal. 
In  view  of  Eq .  (B-10)  ,  Q.  i  s  a  I  so  diagonal;  the  principal  stresses  are 
given  below  in  terms  of  the  principal  strains  and  the  principal  values 
of 


=  X  trer  *  2^iei 

^2  =  ^  *2£2  (B"20) 
r 

°3  =  X  *3  II  0rer  ♦  2/i  ^e3 


The  appropriate  evolution  equation  for  ^  is  obtained  from  Eq.  (B-14) 
and  is  given  in  terms  of  the  principal  components  of  ^  as  follows: 

d*.  =  -  (*.)m 

where 

fk  de .  ,  if  de  -  >  0  ,  e •  >  0 

d*j  =  '  '  .  '  (B-21) 

l  0  ,  otherwise 


We  now  consider  two  special  cases  of  the  above  equations,  namely, 
simple  tension  and  simple  compression. 


[ 

Consider  the  case  of  monoton i ca I ly  increasing  simple  tension 
In  this  case,  <72  =  ^3  =  0,  and  integration  of  Eqs .  (B-21)  with  62  =  £ 
and  the  initial  conditions  $2  =  ^3  =  1/  shows  that  ^2  =  ^3  always 
Hence,  we  find  that 


^2C2  '  ‘  2 (\+fi)  £1^1“  ‘  ^lel  (B-22) 

and 

°i  ■  E*iei 

where 

F  _  (2/i*3\)(i 
11  ~  Ufi 

From  Eq.  (B-21a) ,  the  evolution  equation  for  ^  is 

d^J  +  d^  =  0  ,  (B-25) 


(B-23) 


(B-24) 


subject  to  the  initial  condition  ^ (0) 


1.  Thus,  the  solution  is 


'l 


if  u  s  1 

if  m  M 


(B-26) 


Combining  Eqs.  (B-2.3)  and  (B-26),  23  obtain 


a 


1 


i  f  m  =  1 


,  if  m  *  1 


(B-27) 


Inspection  of  Eq .  (B-27)  reveals  that,  for  0  <  hi  <  1,  ffj  is  not  a 

monotonic  function  of  £^,  but  reaches  a  maximum  and  then  goes  to  zero 
for  finite  £^.  For  1  m  <  3,  0 j  again  reaches  a  maximum  and  then 

decays  monoton i ca I  I y  to  zero  as  fj  4  ».  For  m  =  3,  a j  monoton i ca I  I y 

increases  to  an  asymptotic  constant  value  as  £^  *  ",  while  for  3  <  m  < 
",  increases  monoton  i  ca  I  I  y  with  £^  .  The  influence  of  m  on  the 
constitutive  response  is  shown  in  Figure  B  la.  Thus,  in  the  range  0 
m  <  3,  the  model  exhibits  softening  in  tension  in  agreement  with 
experimental  data. 

The  model  also  unloads  elasticity  with  an  effective  modulus  of 
which  is  diminished  by  the  accumulated  damage  Thus,  the  model 
has  all  of  the  characteristics  of  a  fracturing  elastic  solid.  It 
follows  from  Eq.  (B-22)  that  £2  0  and  hence  $2  =  1>  '•«•,  there  is 

no  damage  in  the  transverse  direction  to  loading. 


Consider  now  the  case  of  monoton i ca I  I y  increasing  uniaxial 
compression.  Here,  0 2  =  0 3  =  0,  while  and  fj  are  compressive 
Again,  Eqs.  (8-21)  together  with  e 2  =  £3,  and  the  initial  conditions 

^2  =  ^3  =  1,  lead  to  the  result  that  ^2  =  ^3  alleys.  Similar  to  the 
case  of  simple  tension,  we  find  that 


*2f2  =  "  ^1£1  (B'28) 

and 

a1  =  (B-29) 


The  damage  evolution  in  this  case  is,  however,  entirely  different  from 
that  for  simple  tension.  Since  both  d£j  and  £}  are  compressive  at  all 
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times,  we  find  from  Eq.  (B-21)  that  d^  =  0  and  thus  ^  =  1  always. 
However,  in  view  of  Eq .  (B-28)  ,  £2  IS  no*  tensile  and  hence  ^2  (an^ 
^3)  will  increase  with  the  consequence  that  damage  will  develop  on 
planes  which  are  parallel  to  the  axis  of  compression,  in  accordance 
with  experimental  observation;  this  is  the  so-called  axial  splitting 
mode  (Horii  and  Nemat-Nasser ,  1985). 


From  Eq.  (B-21),  we  can  write 


d02  ♦  *2k  d€2  =  0 


the  solution  of  which  is 


(B-30) 


e‘ke2 


,  m  =  1 


(B-31) 


(l  -  (l-m)ke2}1“m  ,  m  t  1 


Thus,  in  view  of  Eq.  (B-29)  and  the  condition  ^  =  1,  the  axial  stress 
strain  relation  is 


°i  =  E  ei  ■ 


(B-32) 


so  that  the  response  in  the  axial  direction  is  purely  elastic. 

To  find  the  relation  between  the  strains,  we  use  Eqs .  (B-28)  and 
(B-31)  together  with  the  fact  that  =  1.  It  then  follows  that 


,  m  =  1 


vtl  = 


f 2[l  -  (1-m) k e 2^ 1  m  ,  m  t  1 


(B-33) 
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Here,  we  may  distinguish  three  different  cases: 


(a)  0  ^  m  <  1 

The  strain  is  not  a  monotonic  function  of  £ 2  It  reaches  a 
maximum  absolute  value  and  then  goes  to  zero  for  a  finite  value  of  £2- 
Spec i f i ca I  I y , 

1  (  "\  n+1 

=  b  fe)  <B-34> 

where  n  =  1/ (1-m) . 


(b)  1  <  m  <  2 

Again,  the  strain  £j  is  not  a  monotonic  function  of  £2-  It 
reaches  a  maximum  absolute  value  again  given  by  Eq.  (B-34)  for  1  <  m 
and  then  goes  to  zero  for  infinite  £2-  For  m  =  1 


(B-35) 


where  e  denotes  the  base  of  the  natural  logarithm. 

(c)  2  <  m  <  » 

In  this  case,  £^  is  a  monotonic  function  of  £2  and  approaches 
the  following  limiting  condition  for  increasing  lateral  strain: 


I  im 

£2*®  £1  max 


l  «  ,  m  >  2 


(B-36) 


The  above  cases  are  depicted  graphically  in  Figure  B.l(b). 


Collapse  of  a  Block  Under  Axial  Compression 


It  follows  from  the  above  discussion  that  a  block  consisting  of 
elastic  fracturing  material  will  collapse  under  axial  compression  if  0 
m  2,  since  has  an  upper  bound  given  by  the  expression: 


1 1  max 


J. 

l/k 


2-m 

1-m 


The  limiting  cases  are 


'1 1  max 


'  1 
t/ke 
_1 

.i/k 


m  =  1 
m  =  2 


(B-37) 


(B-38) 


The  collapse  stress  may  be  calculated  directly  from  Eq.  (B-32) ,  with 
the  result: 


2-m 

I  E_  f  1  ] 1-m 

llmax  i/k  12-mJ 


(B- 39) 


Thus,  the  theory  predicts  the  collapse  of  a  block  under  axial  compres¬ 
sion  due  to  damage  on  planes  parallel  to  the  axis  of  compression. 
Again,  this  is  the  axial  splitting  mode. 

Ratio  of  Collapse  Stresses  in  Tension  and  Compression 

At  this  point,  it  is  of  interest  to  explore  the  rationality  of 
the  above  results  by  comparing  the  ratio  oc/o **,  where  oc  and  a ^  are 
collapse  stresses  in  compression  and  tension,  respectively,  obtained 
analytically  with  that  observed  experimentally  for  concrete,  even 


though  the  end-grip  conditions  in  the  tests  may  not  be  ideal,  as 
assumed  in  solutions  obtained  here.  Restricting  attention  to  the  case 
m  =  1,  it  follows  from  Eqs .  (B-37)  and  (B-39)  that 


Since  v  z  0.2  for  concrete,  we  have 


(B-40) 


10 


(B— 41) 


which  is  close  to  experimental  observation  (Raphael,  1984). 
Simple  Shear 


Consider  now  the  case  of  monoton i ca I  I y  increasing  shear  strain 
under  conditions  of  simple  shear  and  small  strain.  The  configuration 
of  interest  is  depicted  in  Figure 

B.2.  Inasmuch  as  the  only  non-zero  strain  component  is  £j2<  we  can 
wr  i  te 

*  2  *12*12  (B-42> 

and 

*  (*il*j2  *  *i2*jl)e12  <B'43) 

Therefore,  in  this  case,  the  general  constitutive  relation 


ffij  =  X  + 


^  *ik*jfLCk l 


(B-44) 


reduces  to  the  form 


a .  . 

'  J 

=  2 

[X  * 

ij^2  +  n 

^ i l^j 2  *  *i 

i  2^j lJ  j 

£  12 

(B-45) 

Therefore, 

,  we 

can 

wr  i  te 

a12 

=  2 

X  j 

12  *  "ki 

^22  +  ^12) ] 

el2 

(B-46) 

all 

=  2[X  ♦ 

^ *11*12 

£12 

(B-47) 

°  22 

=  2(X  ♦ 

2/i]  ^22^12 

C12 

(B-48) 

For 

a  sma 1 1 

stra i n , 

tte  maximum  tensile  strain  will 

lie  in  a 

direction  that  makes  a  45  degree  angle  with  each  of  the  coordinate 
axes.  Consequently,  cracks  will  develop  in  the  material  as  shown  in 
Figure  B.2.  The  components  of  the  unit  normal  n  to  the  crack  are  n^  = 
n2  =  l/JS. 


Since  there  is  only  one  cracking  direction,  we  set  a  =  1  in  the 


evolutionary  equation  for 


'  J 


and  write 


d*ij  =  -  k(<*nrd£  "i"j 


where 


d  =  d.  . n . n . 
”n  r'J  i  J 


de  =  de 


12 


Since  =  n2  =  1/^2*  it  follows  from  Eq.  (B-49)  that 

d^n  =  d^2  =  d^2 


(B-49) 


(B-50) 


(B-51) 


Hence 


(B-53) 


since  ^ • j  =  ^ i j  '"'^iaNy.  From  Eq .  (B-53)  we  find 


^11  =  1  +  ^12  1 

and  Eq.  (B-50a)  gives 


K  =  <*n  +  ** 


12 


(B-54) 


since  according  to  Eq .  (B-51)  . 


22 


Substitution  of  Eq .  (B-53)  into  Eq .  (B-54)  yields 


*  =  1  ♦  2  * 


12 


(B-55) 


which,  when  combined  with  Eq.  (B-49) ,  leads  to  the  result 
M-s  =  -  k(l  .  2*12)" 

Therefore,  since  n^  =  =  1/1*2,  we  can  write 


(B-56) 


d^12  =  ‘  2  +  2  ^12)  d€12 


(B-57) 


which  can  be  integrated  to  give,  for  m  =  1: 
-ke. 


* 


12 


12  . 
e  -  1 


(B-58) 
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and  for  m  t  1 : 

.  \ 

-i- 

¥- * 

ND 

II 

NO  1  >— 1 

1 

k (m-1) € 12+1 

m-1  . 

(B-59) 
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Vi 


V* 

V' 

■  N 


Recalling  Eq.  (B-53)  and  the  fact  that  ^22  =  $ll>  (B-46)  can  be 

placed  in  the  following  form: 


12 


=  2 


(X  ♦  2h)4>12  +  2ji^0  *  n 


12 


12 


(B-60) 


Therefore,  Eqs.  (B-58)  to  (B-60)  describe  the  relationship  between  0^2 
and  £\2' 

From  Eqs.  (B-58)  and  (B-59)  note  that 


“12  *  - 1  <x 


(B-61) 


which  in  view  of  Eq.  (B-60)  leads  to  the  result: 


I  im 
el2‘M” 


a12  =  2  (X  +  2/*)e12 


(B-62) 


Hence,  when  X  <  2ji,  the  ultimate  slope  given  by  Eq.  (B-62)  will  be 
less  than  the  initial  slope,  2 fi,  as  shown  in  Figure  B.3.  The  case  X  > 
2 fi  on  the  other  hand,  leads  to  the  ultimate  slope  being  greater  than 
the  initial  slope,  i.e.,  a  hardening,  which  is  difficult  to  visualize 
from  a  physical  standpoint.  These  limiting  cases,  however,  fall  out¬ 
side  of  the  assumption  of  small  strains  made  here  and  therefore  may 
not  be  physically  meaningful. 


APPENDIX  C 


LIMITING  CONDITION  FOR  HYDROSTATIC  RESPONSE 

INTRODUCTION 

In  Section  4.5,  it  was  pointed  out  that  the  limiting  condition 
at  high  pressures  for  hydrostatic  compression  response  of  concrete  is 


£ 


I  im 
■*e* 


P  P 


a 


00 


(C-l) 


where  a  denotes  the  hydrostatic  pressure  and  is  the  limiting 
volumetric  strain  at  full  compaction.  It  is  shown  in  this  appendix 
that  the  endochronic  model  described  in  Section  4.5  satisfies  this 
condition  when  the  hardening  function  is  taken  in  the  form: 


e* 

P 

e*-e 

P 


PJ 


n 


(C-2) 


and  appropriate  restrictions  are  placed  on  n(>0) .  The  following  proof 
of  this  statement  is  due  to  Valanis  (1986). 

The  Model 


We  consider  the  following  endochronic  model  of  hydrostatic 
response: 


a  = 


de 

K(z  "  z')  dz 

0 


(C-3) 
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where  K(z)  is  weakly  singular  and  integrable, 


dj  =  Ide, 


W 


We  consider  in  the  form: 


p  - c - 

rH  £  *  -  c 
P  P 


,  n  >  0 


and,  in  addition,  consider  the  case  where 
K  (z)  =  ,  0  <  X  <  1 


in  which  $(z)  is  bounded  from  above  and 
I  im  4>(z)  =0 

Z+00 


Ana lysi  s 


We  begin  by  rewriting  Eq.  (C-3)  in  the  form 


[*(?) 

0 


-  *'(<r)]x 


. 


which  under  monotonic  straining  conditions  becomes 


in  view  of  Eq.  (C-4) .  Since  $  is  bounded  from  above  and  is  a  mono 
tonically  decreasing  function  of  its  argument,  we  can  write 


4>(z  -  z)  2  $(z) 


(C-ll 


so  that 


a  >  *(z)  - ^ - r 

[z(f)  -  z'(?')3 
J0 


(C-12 


Thus,  to  show  that  condition  (C-9)  is  satisfied,  it  suffices  to  show 
that 


elxh  [*(<r)  -  z  (<r')]x 


(C-13 


To  show  this,  we  use  Eqs.  (C-4)  and  (C-5)  to  obtain  the  expression: 


C*  -  £ 

dz  =  d? 

P 


which  may  be  integrated  to  give 

z  _  - 1 -  [f£*]n+1  ~  fe*  "  ?]n+1l  > 

(n+1) (e*) n  11  PJ  1  P  J  J 


(C-14 


(C-15 


where  the  condition  z (0)  =  0  has  been  observed.  We  may  therefore 

write  that 


(z-z  ) 


AC 


[(Tor1-  t^r1/ 


Upon  introducing  the  transformation 


X  =  -  £  .  y#  =  e*  -  C' 


we  can  rewrite  Eq.  (C-16)  in  the  form: 

r* 

X ,  ^  nX 


(z  -  z  ) 
0 


£  * 
P 


AC. 


f(y') 


.xn«-l  ^n+ljX 


In  view  of  Eqs .  (C-17) ,  it  then  follows  that 


P  P 


— ^  •  ("*1)x(£5)nX 

(z-z  ) 


e* 

P 


dy' 


(y')x(n+1) 


which  can  be  integrated  to  give 


p  p 


dc' 


(z-z')X 


(n4l)Xfe*'|nX  ,  w  .. 

v  l  p-1  rv'^1_x(n+1) 

[X(n*l)  -  1]  } 


e* 

P 


An  inspection  of  the  righthand  side  of  this  equation  reveals 
becomes  infinite  when 


(C-16) 


(C-17) 


(C— 18) 


(C— 19) 


(C-20) 


that  it 


X (n*l)  -1  <  0 


(C-21) 


or 

n  >  l  -  1  (C-22) 

The  condition  (C-22)  gives  the  relation  between  n  and  X  that  must  be 
satisfied  by  the  hardening  function  (C-6)  in  order  to  satisfy  the 
limiting  condition  (C-l) . 


APPENOIX  D 


DESCRIPTION  OF  CONCRETE  UNDER  COMPRESSIVE 
STRESS  STATES 

1 .  INTRODUCTION 

It  is  well  known  that  the  strength  of  concrete  increases  with 
increasing  mean  confinement  stress.  This  effect  can  be  observed  in 
the  triaxial  failure  surface  of  concrete  in  stress  space,  Figure  D.l. 

In  addition  to  strength,  the  duct i I i ty  of  concrete  is  observed 
to  increase  with  increasing  mean  compressive  stress.  Typical 
increases  in  ductility  during  standard  triaxial  tests  are  shown  in 
Figure  D.2.  The  effect  of  confinement  on  ductility  is  seen  to  be  much 
greater  than  that  on  strength.  This  can  also  be  deduced  from  the 
behavior  of  reinforced  concrete  in  the  presence  of  confinement  steel, 
Figure  D.3. 

Analytical  and  experimental  studies  (see  Nemat-Nasser ,  et  aL., 
[1])  on  brittle  materials  such  as  concrete  suggest  that,  at  low 
confining  stress  levels  failure  occurs  by  the  growth  and  connection  of 
a  narrow  zone  of  microcracks  to  form  one  or  at  most  several 
macrocracks  or  macrofaults,  Figure  D.4a,b.  The  onset  of  this  event 
(formation  of  a  macrocrack)  corresponds  to  the  onset  of  "strain 
softening"  [2].  In  contrast,  at  higher  levels  of  confining  stress 
microcracks  grow  in  a  more  uniform  manner  and  a  distribution  of  micro- 
cracks  evidently  develops,  Figure  D.4c.  Although  failure  may  again 
occur  via  localization,  the  formation  of  this  d i str i but i on  appears  to 
be  responsible  for  the  observed  increase  in  ductility. 


AXIAL  STRESS  (ksl) 


AXIAL  STRAIN 


Figure  D.2.  Influence  of  confinement  on  strength  and  ductility 


jf  reinforced  concrete  in  compress i 


concrete  in  compress 


Below  the  brittle-ductile  transition,  slip  along  the  planes  that 
constitute  the  above  distribution  is  a  more  plausible  explanation  of 
increased  ductility  than  is  true  plastic  flow.  consequently, 
conventional  phenomenological  e I astop I ast i c  descriptions  of  concrete 
under  moderate  compressive  stress  states  appear  to  be  inappropriate. 

In  view  of  the  above  discussion,  an  effort  was  made  to  construct 
a  theory  of  concrete  based  on  a  "slip  system"  concept.  A  description 
of  progress  made  to-date  on  this  subject  follows. 

2.  THE  SLIP  SYSTEM 

Slip  system  models  have  been  used  previously  to  describe  both 
metals  and  frictional  materials.  A  representative  cross-section  of 
such  models  includes  the  works  of  Mandel  [3],  Spencer  [4],  Asaro  [5], 
and  Nemat-Nasser ,  et  aL.,  [6].  An  excellent  discussion  and  in-depth 
treatment  of  the  subject  is  given  by  Nemat-Nasser  [7], 

Figure  D.5  describes  the  slip  concept  in  the  case  of  a  single 
slip  system.  Locally  a  represents  a  unit  normal  to  the  slip  surface 
and  s  is  a  unit  tangent  vector  which  describes  the  d i rect i on  of  slip. 
The  quantity  A  denotes  an  effective  slip  surface  spacing.  The  slip 
system  and  its  geometry  is  presumed  to  form  at  a  certain  critical 
stress  state.  The  collection  of  these  states  forms  a  surface  in 
stress  space. 

In  the  sequel,  a  theory  is  described  for  an  arbitrary  number  of 
slip  systems.  Associated  with  the  ath  system  is  a  normal  vector 
tangent  vector  s^,  and  spacing  A^a^,  Figure  D.5,  where  a  =  1,2  ..., 
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According  to  the 
(i .e. ,  "plastic  si  ip") 
surfaces.  The  material 
elastic. 


slip  concept,  a  I  I  irreversible  deformation 
is  assumed  to  take  place  along  the  slip 
between  these  surfaces  will  be  modeled  as 


3 .  DEFORMATION 


Let  the  relative  rate  of  slip  across  a  slip  surface,  Figure  D  5, 
be  denoted  as  [v^sj].  Further,  let  2;  (i  =  1  to  3)  denote  the  unit 
orthogonal  base  vectors  of  a  rectangular  Cartesian  reference  system. 
Then,  if  the  rate  of  deformation  tensor 


(D-l) 


where  At  is  the  velocity  vector,  is  decomposed  into  elastic  and  inel¬ 
astic  parts  according  to 


d.  . 
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one  obta i ns 


1  (s.n. 

2  l  i  j 


♦  s.n 
J 


where 


s .  =  s  •  e .  n.  =  n  •  e 

I  ~  ~  |  |  **  I 


(D-2) 


(D-3) 


(D-4) 


In  the  case  of  multiple  slip  systems,  the  effects  are  superposed  to 


1 


The  physical  basis  for  (b)  is  the  formation  of  microvoids  and 
microcracks  due  to  slio,  while  that  for  (c)  is  the  collapse  of  voids 
under  hydrostatic  compression.  These  effects  can  be  modeled  by  adding 
to  (D— 7)  terms  of  the  form 


Z  [*?*■«? 


i  j  2  ij  3J  Aa  ij  3 


Consequently,  the  inelastic  deformation  rate  tensor  becomes 


p.  =  jrP*.  4  es. . 

■  i  *  - '  i  i  rt  II 


U  1 J  Aa 


(D-8a) 


where 


S  2  (si"j  *  sj"i)  *  *  Vj92  *  5ij93  '  <D-8b> 


In  addition  to  the  rate  of  deformation  tensor  dfj ,  a  spin  tensor 


u/Pj  is  now  defined  according  to  (Asaro  [5]) 
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where 


_  1  f  a  a  a  a-) 
=  7)  s.n.  -  s.n. 
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(D— 9b) 


4. 


INTERFACE  BEHAVIOR 


The  ath  slip  system  is  assumed  to  be  activated  if  the  resultant 
shear  stress  has  attained  a  critical  value  defined  by 

f a[ra ,oa ,v,6aP)  =  0  ;  a.fi  =  1  to  N  (D-10) 

where  aa  denotes  the  normal  stress  on  the  slip  surface,  p  is  the  mean 
pressure,  and  eaP  are  measures  of  the  history  of  slip.  Slip  occurs 
if,  in  addition  to  (D-10), 

ia[Ta,oa,p,eaP)  =  0  .  (D-ll) 


Equation  (D-ll)  can  be  expressed  in  the  form 


a  La 
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+  aa  tanj?  +  p  tan*79  -  X!  — Mr 

1  2  p=l  v 


(D-12) 


where  tan  rj j  and  tan  rj 2  are  material  parameters  that  characterize, 
respectively,  the  effects  of  normal  stress  and  hydrostatic  stress  on 
slip.  Here  r\ j  corresponds  to  the  usual  friction.  The  quantity  ^2 
reflects  the  influence  of  pressure  on  plastic  flow;  e.g.,  in  crystal 
plasticity  this  term  would  represent  the  influence  of  pressure  on  the 
motion  of  dislocations  £..  (In  two-dimensional  problems  (D-12)  can  be 
re-written  in  terms  of  aa  only,)  The  parameters  kaP  characterize  the 
manner  by  which  the  slip  system’s  resistance  changes  due  to  slip.  In 
what  follows  haP  is  assumed  to  be  symmetric  (but  not  necessarily  posi¬ 
tive  definite).  Finally,  the  effect  of  cohesion  has  been  excluded  in 
the  above  formulation. 


Consider  now  the  problem  of  calculating  the  time  rate  of  change 
of  r®  and  aa  in  (D-10)  -  (D-12) .  Since 


a  "a  a  a  a  a  2a  a  a  a 

t  =  T  •  s  =<r..n.s.  ,  a  =  T  •  n  =  a .  .n.n. 

~  *  ij  i  j  ~  ~  ij  i  j 


(D-13) 


it  is  clear  that  one  must  determine  a®  and  $®  in  order  to  calculate  r® 
and  aa .  The  former  rates  must  be  objective.  Based  on  the  physical 
processes  under  consideration,  the  following  rates  of  change  of  the 
unit  vectors  a®  and  sa  are  defined 


•a  e  a  a  e  a 

n.  =  w. .n.  ,  s.  =  w. -s. 

i  ij  J  •  'J  J 


(D— 14) 


According  to  (D-14) ,  the  unit  vectors  are  corotational  with  the  elas¬ 
tic  deformation  of  the  slip  system.  With  use  of  (D-14),  equations  (D- 
13)  furnish 


a  V  a  a  'a  V  a  a 
t  =  a  ■  .  n .  s  “  ,  o  =  a .  .  n .  n . 
•J  •  J  'J  '  J 


(D— 15) 


where 


$. .  =  a . .  -  we.  a.  -a,  ■  -  w?,  a.  . 
ij  ij  ikkjki  jkki 


(D— 16) 


If  one  defines  the  quantity  q^  according  to 


a.lfaa  a  a)  a  a  .  .  /rv 

q..  =  o  s.n.  ♦  s.n.  ♦  n.n.  tan  n,  0-  ■  tan  (D-17) 

^ij  2  l  i  j  j  iJ  i  j  '1  ij  '2  v  ' 


then  (D-12)  furnishes  the  following  relation  which  can  be  used  to 
ca leu  I  ate  $i j : 


5.  GLOBAL  CONSTITUTIVE  RELATIONS 


In  this  section  the  global  constitutive  relation  for  ajj  is 


assembled.  For  this  purpose  the  elastic  part  of  the  slip  system 
deformation  is  postulated  in  the  form 


l.  = 


L- ..  «d. 


ij  ‘  ijkrkt 


(0—19) 


where  Ljjk£  is  the  instantaneous  elastic  tensor  modulus.  If  the  mate¬ 
rial  between  slip  planes  is  isotropic,  then 


Li jke.  =  G[5ik6j«.  *  ffi^6jk)  *  x<5i j6ki 


(D-20) 


It  wi  I  I  be  assumed  that  Ljjj^  is  independent  of  the  rate  of  elastic 
deformation . 


With  use  of  (0-2),  (D-8a) ,  and  (D— 18) ,  the  relation  (D-19)  can 
be  cast  in  the  following  form: 


^i j  =  C*j k£  (dk£  "  5kt 


(0-21) 
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6 .  REMARKS 

The  constitutive  relations  (D-21)  are  fully  nonlinear  and  incor¬ 
porate  an  arbitrary  number  of  slip  systems.  Their  form  is  identical, 
for  all  practical  purposes,  to  the  relations  obtained  by  Nemat-Nasser 
[7],  The  only  relevant  difference  concerns  the  interpretation  of  the 
tensor  Ljjj^.  The  definition  of  this  quantity,  which  may  depend  on 
the  current  stress  and  the  history  of  deformation,  is  somewhat  vague 
for  brittle  materials  such  as  concrete. 

The  constitutive  relations  (D-21)  include  classical  plasticity 
as  a  special  case.  However,  there  is  an  important  difference.  It  can 
be  easily  shown  that  the  current  model  is  such  that  the  stress  tensor 
a jj  is  noncoaxial  with  the  inelastic  deformation  rate  tensor  d^j . 
This  feature  furnishes  an  ability  to  model  localization  phenomena. 

A  number  of  special  case  studies  using  the  foregoing  results  are 
under  current  study.  These  investigations  are  intended  to  determine 
the  modeling  capabilities  and  limitations  of  the  slip-system  descrip¬ 
tion  of  concrete  behavior. 
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